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We show that a multiply-quantized vortex line core can be energetically stable in a harmonically 
trapped two-component atomic Bose-Einstein condensate. The structure, in which the condensate 
component with the vortex line is surrounded by the second component forming a vortex ring, 
can be identified as a particlelike Skyrmion with a multiply-quantized winding number, and closely 
resembles cosmic vortons. 

PACS numbers: 03.75.Lm,03.75.Mn,12.39.Dc 



Among the most remarkable experiments in atomic 
Bose-Einstein condensates (BECs) are the studies of 
topological defects, such as quantized vortex lines and 
vortex lattices 0. A quantized vortex in a single- 
component BEC exhibits a line singularity of a ID order 
parameter field. However, truly 3D topological objects 
can be afforded by multi-component BECs. It was re- 
cently shown that a localized particlelike soliton, or a 
Skyrmion, can be energetically stable in a trapped two- 
species ^^Rb BEC under realistic experimental conditions 
In this paper we demonstrate that also multiply- 
quantized Skyrmions can be stable in a two-component 
BEC when the interparticle interaction between the two 
species is sufficiently strong. In the stable configurations 
one of the BEC components forms a multiply-quantized 
vortex line and is confined in a toroidal region inside the 
second component forming a vortex ring; see Fig. ^ The 
second component provides the stability of the multiply- 
quantized vortex line core against splitting, and no ex- 
ternal quartic potential is required, as in the case of 
multiply-quantized vortices in a single-component BEC 
0. Moreover, the centrifugal energy, associated with the 
circulation of the vortex line, prevents the shrinking of 
the vortex ring to zero size. As first noted in Ref. 
the structure is closely related to cosmic vortons which 
are superconducting cosmic strings and may have been 
formed in the early Universe phase transitions [J]. Also 
nuclear matter in neutron stars can have a similar de- 
scription I^. The analogies provide an interesting pos- 
sibility for cosmology in laboratory experiments using 
atomic BECs, as performed in superfluid liquid helium 
systems H,I3- 

Multi-component BECs can be prepared in magnetic 
traps by simultaneously confining different atomic species 
in the same trap or in optical dipole traps where the 
spin dynamics is no longer constrained by magnetic fields 
Here we consider two BEC components in per- 
fectly overlapping isotropic magnetic trapping potentials. 
Such a system has been experimentally realized using the 
I T) = |F = 2,m/ = 1) and | [) = |1,-1) hyperfine 
spin states of ^''Rb. In that case the inter- (a^i) and in- 
traspecies {a^^ and a^) interaction strengths are nearly 



equal, with : a^i : a^^ :: 1.03 : 1 : 0.97 9], and the 
order parameter of the two-component BEC has an ap- 
proximate SU(2) symmetry. Since the scattering lengths 
satisfy a^^ > a^^a^, the two species experience dynam- 
ical phase separation and can strongly repel each other 
9]. The interatomic interactions of the two ^^Rb com- 
ponents do not mix the atom population and the atom 
numbers of the both species are separately conserved. 

The properties of the particlelike Skyrmions in a two- 
component atomic BEC have been addressed by several 
authors HHIillillilllllii. An analogous config- 
uration may also be formed in a spin-1 BEC 11]. Other 
theoretical studies of spin defects and textures in multi- 
component BECs include, e.g., the structure of point- 



like defects 
18, 19,(2lil2 




and singular and coreless vortices 
12 3j . Vortices have been experimentally 
studied in two-component BECs ^2^ l26j and a core- 
less vortex has recently been engineered in an optically 
trapped ^^Na spin-1 BEC experiment j^^. Moreover, 
topological defects in a two-component fermionic atomic 
gas in an optical lattice can result in particle number 
fractionilization ^2^. 

A Skyrmion is a localized particlelike soliton such that 
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FIG. 1: Density profiles for energetically stable multiply- 
quantized vortex cores in trapped Skyrmions. We display 
the ID densities of the two BEC components -0- (top) and 
-0+ (bottom) along the x axis (in units of 10~^/~^). We show 
the vortex line core of with the winding number one (solid 
line), two (dashed line), and three (dashed-dotted line) with 
the corresponding density profiles of -0+ exhibiting a singly- 
quantized vortex ring on the xy plane. The size of the vortex 
line core region is approximately 2/, 3/, 4.5/, respectively. 
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the order parameter field becomes uniform, with a con- 
stant asymptotic value, at sufficiently large distances 
from the particle 0|. As a result, we may enclose the 
Skyrmion by a sphere on which the field configuration 
has the constant asymptotic value; with the whole sur- 
face of the sphere identified with a single point. Then the 
physical space inside the sphere can be considered as a 
3-sphere (a sphere in 4D). Moreover, the order param- 
eter of a non-singular two-species BEC takes values on 
[or on SU(2)], and we may assign a winding number 
for the Skyrmion to characterize the mappings from the 
physical space to the order parameter space: ^ S^. 

For the soliton solution we parametrize the or- 
der parameter space by the unit vector n = 
[— Im(?/^_), Re(?/^_), — Im(?/^+), Re(?/^+)]/y^, in terms of 
the real and the imaginary parts of the two condensate 
wave functions '0+(r) and ?/^_(r). It is assumed that the 
total density p{r) = |?/^+(r)p + |?/^_(r)p is nonvanishing 
inside the sphere enclosing the Skyrmion. We then ob- 
tain the nonsingular configuration: 

/?/^+(r)\ _ ^/ cos A — i sin A cos/? \ 
\ijj-{r)J I — i sin A sin/3exp(i7) y ' ^ ^ 

where (A, 7) denote the spherical angles in 4D, such 
that n = (sin A sin cos 7, sin A sin /3 sin 7, sinAcos/3, 
cos A). Moreover, the winding number in terms of the 
spherical angles a = (A, 7) reads: 

W = ^ j d^x sin' Asin/3Det (^|^^ . (2) 

The integration is over the volume of the sphere enclosing 
the soliton, where p is assumed to be nonvanishing, and 
Xj denote the Cartesian coordinates. 

We numerically investigated the energetic stability of 
the BEC wave functions for the initial states which have 
different values of the winding number W . Although gen- 
erally A = A(r) is anisotropic, as an initial state we may 
use the spherically symmetric solution A = A(r), with 
A(0) = 0, since the minimization of the energy does not 
preserve the spherical symmetry of the initial state. It 




FIG. 2: Charge densities for energetically stable trapped 
Skyrmions with W — 1,2,3 (from left to right), for ^ > 0. 
On the xz plane the charge density between isosurface sec- 
tions ranges from W = 0.001 (in units of at the boundary 
(red) to W = 0.017 (left), W = 0.075 (middle), W = 0.055 
(right) at the center (purple or dark grey). The charge is 
concentrated close to the interface of the two species. 



is also illuminating to study the general structure of the 
solutions using this spherically symmetric ansatz. By in- 
serting /3 = kO and 7 = mcj) in Eq. (O, with k and m 
integers and (^, (j)) the spherical angles, we obtain, for 
the even values of /c, W = 0, and for the odd values of k: 

|^=^{A(ro)-isin[2A(ro)]} , (3) 

where ro is the radius of the sphere enclosing the soliton. 
The winding number has an integer value W = nm, for 
'^(^o) = ^TT, with n integer. For A(ro) = (n + l/2)7r, we 
obtain 'hedgehog' solutions with W = {n -\- l/2)m. 

The Skyrmion solution (QJ, with /3 = 6>, 7 = m^, and 
A(0) = and A = nyr at large distances from the particle, 
represents a vortex line, with the winding number m, in 
the component ip- (the "line component") and n vortex 
rings in ip-^ (the "ring component"). The vortex rings 
are formed on the circles z = 0, r = A~^[(n — l/2)7r]. 
The line component, with the vortex line oriented along 
the z axis, is confined inside the ring component and fills 
the toroidal core regions of the vortex rings. The wind- 
ing number W is topologically invariant if the asymptotic 
boundary values are not altered and p remains nonvan- 
ishing for r < tq. We may also obtain a solution for 
W = nm, e.g., by having a single vortex ring in ?/^+, with 
the winding number n, threaded by m singly-quantized 
vortex lines in ?/;_, as long as ip- with m units of circu- 
lation is confined inside the circulating ip-^ . By assuming 
a two-species BEC to be an incompressible gas with a 
constant total atom density in the homogeneous space, 
the stable W = 2 Skyrmions were obtained in Ref. [T^. 
We show that the multiply-quantized vortex cores in a 
two-species BEC tend to split due to the inhomogeneous 
density, similarly to the vortices in a single-component 
BEC 3]. Nevertheless, we demonstrate it is still possible 
to find stable configurations. 

The energetic stability of the Skyrmion in an inhomo- 
geneous finite-sized system is a highly nontrivial concept. 
The configuration may collapse, the asymptotic bound- 
ary conditions of the trapped BEC may be violated, or 
the structure may not remain as singularity free. Gen- 
erally, additional physical mechanisms, such as rotation 
or optical potentials, will be required for stability 0. 
Therefore, it is remarkable that energetically stable par- 
ticlelike solitons can exist under experimental conditions 
of ^^Rb, even in the presence of a strong external con- 
finement due to the harmonic trap 

Even in an infinite, homogeneous system topological 
particlelike solitons can be unstable against shrinking to 
zero size. If the energy does not contain any higher order 
gradient terms than quadratic, a simple spatial rescaling 
of the wave functions shows that the energy decreases 
monotonically with the size of the soliton, in the ab- 
sence of additional stabilizing features. However, in a 
two-component BEC, without internal spin dynamics, 
the number of atoms in the two species is separately con- 
served, which can stabilize the Skyrmion against collapse 
to zero size ITsj. In the phase separation regime the two 
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species can strongly repel each other and the toroidal fill- 
ing due to the line component provides a 1 /r^ centrifugal 
barrier that can prevent the vortex ring from contracting 
to zero radius due to its tension. 

In Ref. 2] it was numerically found that the phase 
separation of the two species is necessary, but not suf- 
ficient, condition for the stability of the Skyrmion with 
the winding number ly = 1 in a trapped BEC. Using the 
parameters of the ^''Rb vortex experiments it was 
also required for the stability that the minimum number 
of atoms is several millions 0- If, in addition, the two- 
component BEC was rotated at the speed in a narrow 
window close to O-OSScj, where uo denotes the trap fre- 
quency, the Skyrmion was stable, without the vortex line 
drifting in the inhomogeneous trap towards the edge of 
the atomic cloud. Alternatively, if the BEC density pro- 
file remained cylindrically symmetric, the stability could 
be provided by rotating sufficiently fast only the vortex 
line component. The stable particlelike soliton studied 
in Ref. [3 consisted of a singly-quantized vortex ring in 
^''Rb I I) component, threaded by a singly-quantized 
vortex line in | |). In this paper we show that also 
multiply-quantized Skyrmions can be energetically sta- 
ble in a trapped BEC, providing a unique system with 
stable multiply-quantized vortex lines. 

We numerically studied the energetic stability of the 
particlelike solitons using the 3D nonlinear mean-field 
theory of the trapped two-component BEC. We mini- 
mized the energy of the Skyrmions by evolving the cou- 
pled Gross-Pitaevskii equations (GPEs), 



ihdtil^i = {Ho -^22^' 
3 



(4) 



in imaginary time, for initial states exhibiting different 
winding numbers W, using the split-step method on a 
spatial grid of 128^. Here Hq = -h'^V'^ /2m^V{r), with 
a perfectly overlapping isotropic trapping potential for 
both states V(r) = mcj^r^/2, and hZij = Airh^aijN /m^ 
where N denotes the total atom number. We assume 
the harmonic oscillator length / = (h/muY^'^ = 3.86/im, 
which for ^^Rb corresponds to :^ 27r x 7.8 Hz j2J, 

®®® 

FIG. 3: The decay of a doubly-quantized vortex line core. 
The density profiles IV^-I^ (for 2/ > 0) show the two singly- 
quantized vortices on the xz plane between isosurface sections 
ranging from low at the boundary (red) to high at the cen- 
ter (purple or dark grey) after the double vortex has already 
split in the imaginary time evolution. On the left, the con- 
stant surface plot of (red torus) represents a vortex 
ring encircling the both vortex lines, and the Skyrmion wind- 
ing number W — 2. In the middle and on the right, at later 
times, one of the vortices has crossed the vortex ring which 
now only encircles one vortex line, and W — 1. 



The winding number was calculated during the imag- 
inary time evolution to determine the stability of the 
Skyrmion. We also numerically minimized the energy 
for the BECs rotating along the z axis, in order to sta- 
bilize the Skyrmion against the drift towards the BEC 
boundary where the density is lower. This was done 
by evolving Eqs. in imaginary time, in the rotat- 
ing frame, obtained by replacing by — VtLz^ with 
Lz = ih{ydx — xdy) and the rotation frequency denoted 
by r^. At every time step, we separately normalized the 
wavefunctions to Ni/N to fix the atom number in each 
component. The simulations were fully 3D, without im- 
posing any symmetry on the solution as it relaxed. 

For particlelike solitons with multiply-quantized vor- 
tex lines in the line component, we used as an initial 
state Eq. Q with W = nm^ for n = 1 and for differ- 
ent integer values of m, embedded in a parabolic den- 
sity profile. As shown before, we obtain this by choosing 
/3 ^ 6>, 7 ^ m6>, and a monotonic function A(r), with 
A(0) = and A = tt at the BEC boundary. In order 
to obtain stable Skyrmion configurations, N had to be 
sufficiently large and the population difference between 
the two species |A^+ — A^_| sufficiently small. In the ex- 
perimental values of ^^Rb the slight difference between 

the scattering lengths a++ : a+_ : a :: 1.03 : 1 : 0.97, 

with a+_ :^ 5.5nm 9], indicates that the SU(2) order 
parameter symmetry is only approximate, being broken 
down to U(l)xU(l). As a result of the phase separa- 
tion for > a++a , the scattering length difference 

was previously found to provide a sufficient repulsion be- 
tween the two species to stabilize the W = 1 Skyrmion, 
when N > 8 x 10^, for = 7V_ and / = 3.86/im 
0, IH^. However, here we were unable to find stable 
multiply-quantized vortex line cores for these parame- 
ters by considering N up to 10^ atoms. For large 
when the Skyrmions did not collapse via shrinking, the 
multiply-quantized vortex cores for m > 1 split into 
singly-quantized vortex lines which drifted towards the 
edge of the BEC. The winding number W decayed by 
integer-multiples always when the vortex lines crossed 
the vortex ring core, by forcing the total density to zero 
at the crossing point and, therefore, by creating a singu- 
larity in the spin texture. In Fig. |31 we show the decay 
of a W = 2 Skyrmion when the doubly-quantized vor- 
tex line has split into two singly-quantized vortices drift- 
ing to opposite directions. One of the vortices crosses 
the vortex ring core, which is numerically also confirmed 
by the sudden change in the value of the winding num- 
ber from two to one. Here TV :^ 10^, = 27V_, and 
Q 0.09/a;. Increasing the rotation speed did not stabi- 
lize the Skyrmions with m > 1, since additional vortices 
entered the ring component due to faster rotation. 

We found energetically stable multiply-quantized vor- 
tex line cores for < 10^ when we considered a stronger 
repulsion between the two BECs with a+_ = a++ o:^ 
5.67nm and a-- 5.34nm. In Fig. ^ we show the stable 
density profiles of vortex lines with the winding numbers 
m = 1,2,3, representing Skyrmions with W = 1,2,3, 
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for the case = N_, N 7.8 x 10^, a+_ 5.89nm, 
a++ 5.67nm, and a__ :^ 5.34nm. Applying the rota- 
tion Q 0.09/cj was found to be sufficient to stabihze 
the configurations at the trap center. The correspond- 
ing topological charge densities W(r), determined from 
Eq. with W = J (i^rW(r), are displayed in Fig. El 

It is also interesting to consider multiply-charged par- 
ticlelike solitons with n > 1 in the case of one singly- 
quantized vortex line (m = 1). We found the stability 
to be very sensitive to the initial conditions. For two 
vortex rings on the xy plane with slightly different radii, 
we found energetically stable W = 2 configurations even 
using the scattering length values of the present ^^Rb ex- 
periments 0- In Fig. El we show the stable case using the 
same parameters as in Fig.Q] However, with one doubly- 
quantized vortex ring as an initial state, the double-ring 
eventually decayed to one singly-quantized ring, with the 
other ring pinching off at the center. This may explain 
why stable W = 2 states with one vortex line were not 
found in the homogeneous space in Ref. where only 
a multiply-quantized vortex ring was considered. 

The particlelike solitons could be created by a sequence 
of Raman pulses, with three orthogonal standing waves 
imprinting a vortex ring, as proposed in Ref. A 
doubly-quantized vortex line has recently been experi- 
mentally prepared at the MIT by changing the magnetic 
field direction. Alternatively, multiply-quantized vortices 
could also be phase-imprinted as proposed in Ref. jsij . 



In conclusion, we showed how multiply-quantized vor- 
tex lines can be energetically stable, forming particlelike 
solitons with multiple winding numbers. The studied set 
of interaction parameter values could possibly be experi- 
mentally reached for different two-species atom mixtures, 
e.g., by using the Feshbach resonances. Moreover, we 
found the stable W = 2 solitons for two vortex rings using 
the parameters of current ^^Rb experiments. Although 
the stability requires reasonably large atom numbers, the 
stable configurations could possibly also be obtained us- 
ing a smaller number of atoms, e.g., by ffist loading atoms 
to a weak trap and then adiabatically increasing the trap- 
ping frequencies to obtain the required nonlinearity 30]. 
We did not study the stability of truly giant vortex line 
cores with m ^ 1. This would be interesting, as giant 
vortex systems have also been created experimentally Q . 

This research was supported by the EPSRC. 
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FIG. 4: Density profiles for energetically stable W = 2 parti- 
clelike solitons with two vortex rings. We show the ID den- 
sities along the x axis (in units of 10~^/~^) for ip- with a 
vortex line at x = (dashed blue line), ip-^ (dotted red line) 
with two vortex rings, and the total density (solid line). 
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